This week’s data analysis mostly deals with the human-dataset. The data contained a number of variables mapping human development and gender equality in 195 countries. This data was then cut down to eight variables. They are the following variables measuring development and gender equality in 195 countries:
+ Life.Expectancy refers to a citizens average life expectancy at birth
+ Education.Years is the amount of years a citizen is planned to spend in education
+ GNI.Per.Capita is the Gross National INcome adjusted for population
+ Maternal.Mortality is a ratio of deaths per 100,000 births
+ Teen.Birth is the number of births per 1,000 women ages 15 to 19.
+ Parliamentary.Participation is the ratio of women to men in parliament
+ Female.Secondary.Education is the percentage of females that attend secondary education
+ Female.Work is the percentage of females that participate in the labour force
Below are simple bar-plots graphically showcasing the variables used this week. As the graphs show, most of these variables do not follow a gaussian distribution, meaning that they are not bell-curves. The only exception to this is the variable mapping the planned years a citizen is supposed to spend in school.
Even further down, the reader can find the correlation plot, where "*" indicates statistical significance to the 0.05-level. As can be seen, three sets of variables emerge from the data:
+ Female.Work and Parliamentary.Participation, which are not strongly correlated with anything
+ Maternal.Mortality and Teen.Birth, which are strongly positively correlated only with one another.
+ The remaining variables on education, GNI, and life-expectancy, which are (semi) strongly positively correlated with each other and strongly negatively correlated with the preceding set of variables.
## Female.Secondary.Education Female.Work Life.Expectancy Education.Years
## Min. : 0.90 Min. :13.50 Min. :49.00 Min. : 5.40
## 1st Qu.: 27.15 1st Qu.:44.45 1st Qu.:66.30 1st Qu.:11.25
## Median : 56.60 Median :53.50 Median :74.20 Median :13.50
## Mean : 55.37 Mean :52.52 Mean :71.65 Mean :13.18
## 3rd Qu.: 85.15 3rd Qu.:61.90 3rd Qu.:77.25 3rd Qu.:15.20
## Max. :100.00 Max. :88.10 Max. :83.50 Max. :20.20
## GNI.per.Capita Maternal.Mortality Teen.Birth
## Min. : 581 Min. : 1.0 Min. : 0.60
## 1st Qu.: 4198 1st Qu.: 11.5 1st Qu.: 12.65
## Median : 12040 Median : 49.0 Median : 33.60
## Mean : 17628 Mean : 149.1 Mean : 47.16
## 3rd Qu.: 24512 3rd Qu.: 190.0 3rd Qu.: 71.95
## Max. :123124 Max. :1100.0 Max. :204.80
## Parliamentary.Participation
## Min. : 0.00
## 1st Qu.:12.40
## Median :19.30
## Mean :20.91
## 3rd Qu.:27.95
## Max. :57.50
Below, Principal Component Analysis, relying on Singular Value Decomposition, is applied to the non-standardized “Human”-data explored above. This analysis is expected to turn out somewhat skewed, as PCA relies on variance to identify the Principal Components. Consequently, variables with larger average values will have an unreasonably high impact on the analysis.
## Standard deviations (1, .., p=8):
## [1] 18544.172057 186.283543 25.972416 20.074641 14.321634
## [6] 10.634338 3.721257 1.428190
##
## Rotation (n x k) = (8 x 8):
## PC1 PC2 PC3
## Female.Secondary.Education -9.317458e-04 0.085169068 -0.3652346210
## Female.Work 9.124960e-05 -0.031442122 -0.0180958993
## Life.Expectancy -2.815825e-04 0.028224407 -0.0170864406
## Education.Years -9.562916e-05 0.007556395 -0.0210421918
## GNI.per.Capita -9.999828e-01 -0.005830843 0.0006412388
## Maternal.Mortality 5.655739e-03 -0.987500960 -0.1481355205
## Teen.Birth 1.233962e-03 -0.125305410 0.9184673154
## Parliamentary.Participation -5.526462e-05 0.003211742 -0.0038388487
## PC4 PC5 PC6
## Female.Secondary.Education -0.8435797499 -3.770012e-01 -6.083913e-02
## Female.Work -0.3928157700 8.233860e-01 4.077784e-01
## Life.Expectancy -0.0212996883 1.136966e-02 -6.669649e-02
## Education.Years -0.0308785262 3.896982e-03 -2.437473e-02
## GNI.per.Capita 0.0002381021 6.651486e-06 2.062449e-05
## Maternal.Mortality -0.0173448186 -3.955532e-02 -2.010447e-02
## Teen.Birth -0.3475453954 -1.381061e-01 -2.499459e-02
## Parliamentary.Participation -0.1075784134 3.989026e-01 -9.077136e-01
## PC7 PC8
## Female.Secondary.Education 0.0303039622 3.118655e-02
## Female.Work -0.0093667016 6.654689e-03
## Life.Expectancy -0.9901494274 1.161211e-01
## Education.Years -0.1134676761 -9.925031e-01
## GNI.per.Capita 0.0001028639 1.871381e-05
## Maternal.Mortality -0.0244043639 -7.101321e-04
## Teen.Birth -0.0127951595 -8.079546e-03
## Parliamentary.Participation 0.0704544742 1.925677e-02
## Importance of components:
## PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8
## Standard deviation 1.854e+04 186.2835 25.97 20.07 14.32 10.63 3.721 1.428
## Proportion of Variance 9.999e-01 0.0001 0.00 0.00 0.00 0.00 0.000 0.000
## Cumulative Proportion 9.999e-01 1.0000 1.00 1.00 1.00 1.00 1.000 1.000
As expected, the graph can tell us very little, and PC1 purports to explain (essentially) all of the variation. Let us standardize the data and try again.
By standardizing the variables in the “Human”-data, the PCA is able to work on variables with comparable variances. Accordingly, variables on scales where high values are expected will not have a disproportionately high impact on the calculations. As can be seen below, the biplot produced by the standardized data can offer valuable insights into the data.
## Standard deviations (1, .., p=8):
## [1] 2.1359720 1.1170602 0.8767066 0.7202736 0.5530904 0.5299143 0.4493522
## [8] 0.3372776
##
## Rotation (n x k) = (8 x 8):
## PC1 PC2 PC3 PC4
## Female.Secondary.Education -0.39874523 0.101737848 -0.210481695 -0.31033900
## Female.Work 0.14439668 0.682158952 -0.575062704 -0.27300332
## Life.Expectancy -0.43080653 0.003997077 0.073302641 -0.02783772
## Education.Years -0.41858363 0.138149663 -0.073869337 -0.07719277
## GNI.per.Capita -0.33914119 0.098797891 -0.338769060 0.84020987
## Maternal.Mortality 0.41991628 0.123208094 -0.145471957 0.28520785
## Teen.Birth 0.40271826 0.088902996 -0.003213286 0.11830579
## Parliamentary.Participation -0.07626861 0.687286162 0.691544283 0.14537131
## PC5 PC6 PC7 PC8
## Female.Secondary.Education 0.43437742 -0.50478735 0.48033587 0.12578655
## Female.Work -0.22413576 0.23657515 0.01076780 -0.04754207
## Life.Expectancy 0.04920168 0.57970641 0.04122211 0.68415054
## Education.Years 0.30831448 -0.11017948 -0.81713841 -0.13919229
## GNI.per.Capita -0.01249472 0.01564322 0.16310711 -0.16583233
## Maternal.Mortality 0.05493343 -0.43780233 -0.24286926 0.67254029
## Teen.Birth 0.81248359 0.36928738 0.07464930 -0.11761127
## Parliamentary.Participation -0.01724555 -0.11284562 0.09265604 -0.02894539
## Importance of components:
## PC1 PC2 PC3 PC4 PC5 PC6 PC7
## Standard deviation 2.1360 1.1171 0.87671 0.72027 0.55309 0.5299 0.44935
## Proportion of Variance 0.5703 0.1560 0.09608 0.06485 0.03824 0.0351 0.02524
## Cumulative Proportion 0.5703 0.7263 0.82235 0.88720 0.92544 0.9605 0.98578
## PC8
## Standard deviation 0.33728
## Proportion of Variance 0.01422
## Cumulative Proportion 1.00000
As the caption suggests, the data seems to indicate that while participation of women in the labour-force and political decision-making are important components of gender equality, they do not predict the overall (physical) well-being of women. More of this below. The reason why our standardized data is able to provide more insights into this, is due to the standardized means and variations that it contains. PCA creates its principal components (PC) by calculating a line of best fit through the observations of the variables that account for most of the variance (while maintaining the sum-of-squares for the resulting coefficients at =1). The second PC is calculated from the remaining variables similarly, but in such a way that it does not correlate with the first PC. With standardized data, the variations are not affected by the scales of the units of the variables and are as such comparable - no single variable can account for most of the variance due to the large values it takes. Conversely, with the non-standardized data, the variable measuring the per-capita GNI of the countries had the largest explanatory power, since the values that its variation takes are in the tens-of-thousands. The rest of the variables take much lower values (The second highest absolute values are taken by maternal mortality and even that has a range of 0 to 1000), and consequently they seemingly account for less of the variation in the data. This is why we see PC1 correlate almost 100% with the variable GNI.per.Capita in the biplot run on non-standardized variables and why the said variable seems to have a disproportionately high variation among the dataset.
As indicated above, high labour and political participation of women are not short hand for high level of well being among women. We can see that the data is principally divided between countries where there are high levels death at childbirth and adolescent births, and countries where there are high levels of life-expectancy, education and GNI. We might say that this pricipal component mostly maps reproductive health. Indeed, the former two variables are highly correlated between each other, which is not a surprise: giving birth at a young (teen) age creates a greater risk of death at childbirth due to the fact that the mother’s body is not yet fully developed. The latter four variables on the other hand predict lower levels of maternal mortality and adolescent birth, partly because higher levels of education are known to predict a lower fertility rates, and a delayed start to getting children among women, partly because higher GNI predicts better health care, both contributing to a higher life-expectancy.
On the other hand, labour and political participation presuppose neither of these sets of variables - they are almost completely uncorrelated with them. This indicates that women participate in the labour force and parliament even in poorer countries with more traditional conceptions of child-rearing a worser health-care. On the other hand, it perhaps also indicates that most countries perform suboptimally, when it comes to the participation of women in politics and the labour force. As such those two variables are bad predictors of development either in terms of money, education or health.
The tea-dataset, according to the information provided by ?tea, includes three sets of variables: the first 18 relate to tea drinking habits, the next 5 map personal details (with age featuring twice, once as integers and once as a categorical variable with five categories), and the final 13 questions mapping product perception. For this section I have chosen to divide the variables into two sets: the active, habit-variables being fed into the MCA and the supplementary personal detail variables, which will be added as colors into the plot to see if we can find links between these supplementary variables and the two dimensions of variability showcased in the plot. The product-perception variables will be excluded in this analysis. Furthermore, only the categorical age variable will be included.
Below, simple bar graphs have been created alongside the structure and dimensions of the remaining data. The reader is welcome to explor the graphs, but since their interpretation is not a requested part of this analysis, a few short comments should suffice: - We can see that there is a clear distinction between places where tea is drunk (home, friends) and where and when it is not (pub, restaurant, dinner, lunch, work) - People seem to drink Earl Grey frequently - Slightly more women answered the survey.
## [1] 300 36
## 'data.frame': 300 obs. of 36 variables:
## $ breakfast : Factor w/ 2 levels "breakfast","Not.breakfast": 1 1 2 2 1 2 1 2 1 1 ...
## $ tea.time : Factor w/ 2 levels "Not.tea time",..: 1 1 2 1 1 1 2 2 2 1 ...
## $ evening : Factor w/ 2 levels "evening","Not.evening": 2 2 1 2 1 2 2 1 2 1 ...
## $ lunch : Factor w/ 2 levels "lunch","Not.lunch": 2 2 2 2 2 2 2 2 2 2 ...
## $ dinner : Factor w/ 2 levels "dinner","Not.dinner": 2 2 1 1 2 1 2 2 2 2 ...
## $ always : Factor w/ 2 levels "always","Not.always": 2 2 2 2 1 2 2 2 2 2 ...
## $ home : Factor w/ 2 levels "home","Not.home": 1 1 1 1 1 1 1 1 1 1 ...
## $ work : Factor w/ 2 levels "Not.work","work": 1 1 2 1 1 1 1 1 1 1 ...
## $ tearoom : Factor w/ 2 levels "Not.tearoom",..: 1 1 1 1 1 1 1 1 1 2 ...
## $ friends : Factor w/ 2 levels "friends","Not.friends": 2 2 1 2 2 2 1 2 2 2 ...
## $ resto : Factor w/ 2 levels "Not.resto","resto": 1 1 2 1 1 1 1 1 1 1 ...
## $ pub : Factor w/ 2 levels "Not.pub","pub": 1 1 1 1 1 1 1 1 1 1 ...
## $ Tea : Factor w/ 3 levels "black","Earl Grey",..: 1 1 2 2 2 2 2 1 2 1 ...
## $ How : Factor w/ 4 levels "alone","lemon",..: 1 3 1 1 1 1 1 3 3 1 ...
## $ sugar : Factor w/ 2 levels "No.sugar","sugar": 2 1 1 2 1 1 1 1 1 1 ...
## $ how : Factor w/ 3 levels "tea bag","tea bag+unpackaged",..: 1 1 1 1 1 1 1 1 2 2 ...
## $ where : Factor w/ 3 levels "chain store",..: 1 1 1 1 1 1 1 1 2 2 ...
## $ price : Factor w/ 6 levels "p_branded","p_cheap",..: 4 6 6 6 6 3 6 6 5 5 ...
## $ age : int 39 45 47 23 48 21 37 36 40 37 ...
## $ sex : Factor w/ 2 levels "F","M": 2 1 1 2 2 2 2 1 2 2 ...
## $ SPC : Factor w/ 7 levels "employee","middle",..: 2 2 4 6 1 6 5 2 5 5 ...
## $ Sport : Factor w/ 2 levels "Not.sportsman",..: 2 2 2 1 2 2 2 2 2 1 ...
## $ age_Q : Factor w/ 5 levels "15-24","25-34",..: 3 4 4 1 4 1 3 3 3 3 ...
## $ frequency : Factor w/ 4 levels "1/day","1 to 2/week",..: 1 1 3 1 3 1 4 2 3 3 ...
## $ escape.exoticism: Factor w/ 2 levels "escape-exoticism",..: 2 1 2 1 1 2 2 2 2 2 ...
## $ spirituality : Factor w/ 2 levels "Not.spirituality",..: 1 1 1 2 2 1 1 1 1 1 ...
## $ healthy : Factor w/ 2 levels "healthy","Not.healthy": 1 1 1 1 2 1 1 1 2 1 ...
## $ diuretic : Factor w/ 2 levels "diuretic","Not.diuretic": 2 1 1 2 1 2 2 2 2 1 ...
## $ friendliness : Factor w/ 2 levels "friendliness",..: 2 2 1 2 1 2 2 1 2 1 ...
## $ iron.absorption : Factor w/ 2 levels "iron absorption",..: 2 2 2 2 2 2 2 2 2 2 ...
## $ feminine : Factor w/ 2 levels "feminine","Not.feminine": 2 2 2 2 2 2 2 1 2 2 ...
## $ sophisticated : Factor w/ 2 levels "Not.sophisticated",..: 1 1 1 2 1 1 1 2 2 1 ...
## $ slimming : Factor w/ 2 levels "No.slimming",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ exciting : Factor w/ 2 levels "exciting","No.exciting": 2 1 2 2 2 2 2 2 2 2 ...
## $ relaxing : Factor w/ 2 levels "No.relaxing",..: 1 1 2 2 2 2 2 2 2 2 ...
## $ effect.on.health: Factor w/ 2 levels "effect on health",..: 2 2 2 2 2 2 2 2 2 2 ...
Moving on to the actual Multiple Correspondence Analysis. Below, the reader can first find the variable biplot of the MCA, where variables in red denote active variables and variables in green supplementary variables. It is first worth noting that the two dimensions identified by the MCA account together for only 16.7% of the principal inertia, meaning that most of the variation in the data is left unaccounted for. In fact, the summary-statistic shows that increasing the amount of dimensions has little effect on this, with six dimensions only accounting for 36.8% of the principal inertia.
Be that as it may, we will continue to work with these results. The theory that I offer based on the below information is that dimension 1 maps what could be generalised as “tea consumption.” We see the frequency of tea-drinking increase as we move from left to right. Additionally, the variables on the right-hand side of the plot indicate higher levels of tea consumption - to the point that tearooms are visited. The left-hand side showcase answers indicating little interest in tea: the tea is cheap, from the chain store and the people do not generally choose tea as their drink when in restaurants, at home or at friends. Dimension 2 on the other hand maps “tea enjoyment”, or something similar. We see the higher values correspond with upscale tea bought from a teashop in loose-form. The lower variables reflect a more basic approach to tea. But as said, even if correct, these dimensions and their interpretations do not explain a lot of the variance in the data.
Consequently, it is not surprising that when we apply color to the plots to map the respondents by their belonging to different categories in the four supplementary variables, we have a hard time identifying any significant linkages. Age does not seem to be linked to either of the dimensions identified. With sports, there seems to be indication of individuals that do not take part in sport being also lesser tea drinkers, but the indications for this is weak. Similarly, students and elderly seem to drink more tea when it comes to employment status. Perhaps the most clearest linkage is between gender and dimension 1, where females are more prominent at the right-hand side of the scale, perhaps indicating that women tend to be bigger casual tea drinkers, with men more often opting for other drinks, or drinking finer teas, as indicated by the higher amount of men at the top of dimension two.
Overall, these plots are not helpful. Primarily due to their low ability to account for variation in the data, and also, perhaps consequently, due to the absence of other linkages that they might have identified.
##
## Call:
## MCA(X = Tea_Who_And_How, quanti.sup = NULL, quali.sup = 19:22,
## graph = FALSE)
##
##
## Eigenvalues
## Dim.1 Dim.2 Dim.3 Dim.4 Dim.5 Dim.6 Dim.7
## Variance 0.148 0.116 0.094 0.078 0.073 0.072 0.067
## % of var. 9.392 7.335 5.972 4.957 4.635 4.534 4.233
## Cumulative % of var. 9.392 16.727 22.699 27.655 32.291 36.824 41.057
## Dim.8 Dim.9 Dim.10 Dim.11 Dim.12 Dim.13 Dim.14
## Variance 0.065 0.060 0.059 0.058 0.055 0.052 0.050
## % of var. 4.106 3.829 3.719 3.672 3.500 3.319 3.193
## Cumulative % of var. 45.163 48.992 52.711 56.383 59.883 63.202 66.395
## Dim.15 Dim.16 Dim.17 Dim.18 Dim.19 Dim.20 Dim.21
## Variance 0.049 0.046 0.045 0.044 0.039 0.037 0.035
## % of var. 3.094 2.932 2.841 2.785 2.480 2.330 2.210
## Cumulative % of var. 69.488 72.421 75.262 78.047 80.527 82.857 85.067
## Dim.22 Dim.23 Dim.24 Dim.25 Dim.26 Dim.27 Dim.28
## Variance 0.034 0.034 0.029 0.028 0.027 0.026 0.024
## % of var. 2.149 2.132 1.860 1.755 1.690 1.637 1.529
## Cumulative % of var. 87.216 89.348 91.208 92.963 94.653 96.290 97.820
## Dim.29 Dim.30
## Variance 0.019 0.016
## % of var. 1.196 0.984
## Cumulative % of var. 99.016 100.000
##
## Individuals (the 10 first)
## Dim.1 ctr cos2 Dim.2 ctr cos2 Dim.3 ctr
## 1 | -0.573 0.739 0.160 | -0.182 0.095 0.016 | -0.348 0.427
## 2 | -0.381 0.326 0.139 | -0.127 0.047 0.016 | -0.604 1.288
## 3 | 0.095 0.020 0.006 | -0.143 0.059 0.013 | 0.303 0.325
## 4 | -0.634 0.904 0.281 | 0.009 0.000 0.000 | 0.141 0.070
## 5 | -0.138 0.043 0.023 | -0.100 0.029 0.012 | -0.105 0.039
## 6 | -0.734 1.212 0.265 | 0.018 0.001 0.000 | -0.027 0.003
## 7 | -0.084 0.016 0.008 | -0.158 0.072 0.027 | -0.081 0.023
## 8 | -0.262 0.154 0.054 | -0.022 0.001 0.000 | -0.032 0.004
## 9 | 0.198 0.088 0.033 | 0.179 0.092 0.027 | -0.593 1.244
## 10 | 0.344 0.265 0.081 | 0.436 0.546 0.130 | -0.461 0.752
## cos2
## 1 0.059 |
## 2 0.350 |
## 3 0.058 |
## 4 0.014 |
## 5 0.013 |
## 6 0.000 |
## 7 0.007 |
## 8 0.001 |
## 9 0.295 |
## 10 0.146 |
##
## Categories (the 10 first)
## Dim.1 ctr cos2 v.test Dim.2 ctr cos2
## breakfast | 0.231 0.908 0.049 3.835 | -0.221 1.062 0.045
## Not.breakfast | -0.213 0.838 0.049 -3.835 | 0.204 0.980 0.045
## Not.tea time | -0.483 3.622 0.181 -7.360 | 0.116 0.269 0.011
## tea time | 0.375 2.808 0.181 7.360 | -0.090 0.209 0.011
## evening | 0.331 1.332 0.057 4.133 | -0.013 0.003 0.000
## Not.evening | -0.173 0.696 0.057 -4.133 | 0.007 0.001 0.000
## lunch | 0.714 2.656 0.088 5.121 | -0.352 0.824 0.021
## Not.lunch | -0.123 0.457 0.088 -5.121 | 0.060 0.142 0.021
## dinner | -0.795 1.568 0.048 -3.769 | 0.827 2.174 0.051
## Not.dinner | 0.060 0.118 0.048 3.769 | -0.062 0.164 0.051
## v.test Dim.3 ctr cos2 v.test
## breakfast -3.665 | -0.621 10.329 0.356 -10.315 |
## Not.breakfast 3.665 | 0.573 9.535 0.356 10.315 |
## Not.tea time 1.773 | 0.158 0.607 0.019 2.403 |
## tea time -1.773 | -0.122 0.471 0.019 -2.403 |
## evening -0.169 | 0.510 4.994 0.136 6.383 |
## Not.evening 0.169 | -0.267 2.611 0.136 -6.383 |
## lunch -2.520 | 0.356 1.037 0.022 2.552 |
## Not.lunch 2.520 | -0.061 0.178 0.022 -2.552 |
## dinner 3.922 | 0.637 1.586 0.031 3.023 |
## Not.dinner -3.922 | -0.048 0.119 0.031 -3.023 |
##
## Categorical variables (eta2)
## Dim.1 Dim.2 Dim.3
## breakfast | 0.049 0.045 0.356 |
## tea.time | 0.181 0.011 0.019 |
## evening | 0.057 0.000 0.136 |
## lunch | 0.088 0.021 0.022 |
## dinner | 0.048 0.051 0.031 |
## always | 0.060 0.003 0.053 |
## home | 0.013 0.000 0.131 |
## work | 0.099 0.049 0.002 |
## tearoom | 0.344 0.017 0.000 |
## friends | 0.208 0.011 0.160 |
##
## Supplementary categories (the 10 first)
## Dim.1 cos2 v.test Dim.2 cos2 v.test Dim.3 cos2
## F | 0.168 0.041 3.511 | -0.114 0.019 -2.375 | -0.045 0.003
## M | -0.245 0.041 -3.511 | 0.166 0.019 2.375 | 0.066 0.003
## employee | -0.153 0.006 -1.305 | -0.145 0.005 -1.243 | 0.064 0.001
## middle | -0.010 0.000 -0.071 | 0.312 0.015 2.114 | -0.330 0.017
## non-worker | 0.005 0.000 0.045 | 0.175 0.008 1.580 | -0.254 0.017
## other worker | -0.030 0.000 -0.140 | 0.009 0.000 0.041 | 0.063 0.000
## senior | 0.370 0.018 2.325 | 0.047 0.000 0.297 | -0.106 0.001
## student | 0.031 0.000 0.299 | -0.279 0.024 -2.664 | 0.391 0.047
## workman | -0.454 0.009 -1.601 | 0.216 0.002 0.764 | 0.061 0.000
## Not.sportsman | -0.031 0.001 -0.444 | 0.016 0.000 0.222 | -0.036 0.001
## v.test
## F -0.942 |
## M 0.942 |
## employee 0.545 |
## middle -2.238 |
## non-worker -2.285 |
## other worker 0.292 |
## senior -0.667 |
## student 3.732 |
## workman 0.217 |
## Not.sportsman -0.518 |
##
## Supplementary categorical variables (eta2)
## Dim.1 Dim.2 Dim.3
## sex | 0.041 0.019 0.003 |
## SPC | 0.029 0.044 0.066 |
## Sport | 0.001 0.000 0.001 |
## age_Q | 0.008 0.066 0.103 |